PT-symmetry breaking and laser-absorber modes in optical scattering systems 
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Using a scattering matrix formalism, we derive the general scattering properties of optical struc- 
tures that are symmetric under a combination of parity and time-reversal (PT) ■ We demonstrate the 
existence of a transition beween PT-symmetric scattering eigenstates, which are norm-preserving, 
and symmetry-broken pairs of eigenstates exhibiting net amplification and loss. The system pro- 
posed by Longhi 9 , which can act simultaneously as a laser and coherent perfect absorber, occurs 
at discrete points in the broken symmetry phase, when a pole and zero of the S-matrix coincide. 
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Schrodinger equations that violate time-reversal sym- 
metry due to a non-Hermitian potential, but retain com- 
bined PT (parity-time) symmetry, have been extensively 
studied since the work of Bender et. al. [TJ who 
showed that such systems can exhibit real energy eigen- 
values, suggesting a possible generalization of quantum 
mechanics. Moreover, PT-symmetric systems can dis- 
play a spontaneous breaking of PT-symmetry, at which 
the reality of the eigenvalues is lost [TJ [3] . Although PT- 
symmetric quantum mechanics remains speculative as a 
fundamental theory, the idea has been fruitfully extended 
to wave optics [5HS]. The classical electrodynamics of 
a medium with loss or gain breaks T symmetry in the 
mathematical sense, although the underlying quantum 
electrodynamics is of course T-symmetric. PT symme- 
try is maintained in optical systems by means of balanced 
gain and loss regions that transform into one another un- 
der parity; thus the combined PT operation, which also 
interchanges loss and gain, leaves the system invariant. 

We show in this Letter that the scattering behavior of 
a general PT-symmetric system can exhibit one or mul- 
tiple spontaneous symmetry-breaking transitions. This 
result applies to arbitrarily complex PT-symmetric scat- 
tering geometries, whereas earlier works on optical PT 
symmetry breaking were inherently restricted to waveg- 
uide (paraxial) geometries [3H6] for which resonances in 
the propagation direction play no role. In addition, we 
elucidate the properties of certain singular solutions oc- 
curring in such systems, recently studied for special cases 
by several authors [7H9] , where a pole and a zero of 
the scattering matrix (S-matrix) coincide at a specific 
real frequency. A real-frequency pole corresponds to the 
threshold for laser action, while a real-frequency zero im- 
plies the reverse process to lasing, in which a particu- 
lar incoming mode is perfectly absorbed. A device ex- 
hibiting the latter phenomenon, which does not require 
PT-symmetry, has been termed a "coherent perfect ab- 
sorber" (CPA) [10 . A PT-symmetric scatterer, at these 
singular points, can function simultaneously as a CPA 
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and a laser at threshold, as noted by Longhi [S]. The 
present work establishes the CPA-laser points as spe- 
cial solutions in a wider "phase" of PT-broken scatter- 
ing eigenstates. We identify signatures of both the PT- 
breaking transition and the CPA-lasing points, for several 
exemplary and experimentally-feasible geometries. 

S-matrix properties - Consider an optical cavity cou- 
pled to a discrete set of scattering channels, denoted 
fi = 1, 2, • • • . Incoming fields enter via the input chan- 
nels, interact in the cavity, and exit via the output chan- 
nels. For simplicity, we focus on the scalar wave equa- 
tion, which directly describes one-dimensional (ID) and 
two-dimensional (2D) systems. A steady-state scattering 
solution for the electric field, E(r), obeys 

[V 2 +n 2 {f){u J 2 /c 2 )]E(f)=0. (1) 

For amplifying or dissipative media, n(r) is complex. The 
frequency u> is real for physical processes but can be use- 
fully continued to the complex plane. Henceforth, we set 
c = 1. Outside the cavity, E has the form 

m = E [</V< (r, w) + ¥V< ut <f> w )] • ( 2 ) 

Here itlf (r, uj) and it° ut (r, uj) denote the input and output 
channel modes, whose exact form depends on the scat- 
tering geometry (e.g. plane waves, spherical/cylindrical 
waves, or waveguide modes). The complex input and out- 
put amplitudes, ip^ and c/? M , are related by the S-matrix: 

ES/wMVv = ¥V ( 3 ) 

V 

For all values of w, S(oj) is a (complex) symmetric matrix 
[11] : if n(r) and uj are real, it is also unitary. 

We now review the properties of S{uj) for a complex, 
PT-symmetric system [8]. (Henceforth, P will refer to 
any linear symmetry operation such that P 2 = 1, includ- 
ing not only parity, but also ir rotations and inversion.) 
First, observe that the T-operator maps incoming chan- 
nel modes to outgoing ones: 

Tuf(r,uj)=u°^(r,uj*), (4) 
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where T is the anti-linear complex conjugation operator. 
The generalized parity operator P is a linear operator 
acting on scalar fields satisfying 

P </° ut (r, u)=Yl ^X n/ ° ut (f, ") . (5) 

where the system-dependent matrix V mixes the channel 
functions, but never transforms between incoming and 
outgoing channels. Note that (VT) 2 = 1. 

If the system is PT-symmetric, for any solution ([2]) 
there exists a valid solution (PT)E(r) at frequency uj*: 

E [VTvh <tf w *) + CPT1>)v u T(r> • (6) 

Comparing this to ([2| and pi), we conclude that 

(VT)S{u*){VT) = S- X {u). (7) 

This is the fundamental relation obeyed by VT- 
symmetric S-matrices, and we will now show that it has 
important implications for the eigenvalue spectrum. 

Multiplying both sides of ¥f\ by an eigenvector of 
S(lj) with eigenvalue s n gives 

S(lj*) (VTi>„) = \ (VTi> n ) ■ (8) 

S n 

Hence, the inverse of the complex conjugate of any eigen- 
value of S(lu) is an eigenvalue of S(ui*). For real uj, this 
implies that | det S(uj)\ = 1, just as for S-matrices having 
pure T symmetry, which are unitary. Unitarity imposes 
a stronger constraint: each eigenvalue is unimodular, so 
unitary S-matrices do not have poles or zeros for real uj. 
Both T and VT symmetry imply that poles and zeros 
occur in complex conjugate pairs. 

Symmetry-breaking transition - Eq. |7]) is a weaker con- 
straint than unitarity, and can be satisfied in two ways: 
either each eigenvalue is itself unimodular, or the eigen- 
values form pairs with reciprocal moduli. These two pos- 
sibilities correspond to symmetric and symmetry-broken 
scattering behavior. In ID, there are just two scattering 
eigenvectors, so the entire S-matrix is either in the sym- 
metric or broken-symmetry "phase" . (In higher dimen- 
sions, as we will see, the transition can occur in different 
eigenspaces of the S-matrix, so that the system can be 
in a mixed phase; initially we focus on ID.) Let us de- 
note the eigenvectors by ip±. In the symmetric phase, 
each ip± is itself PT-symmetric, i.e. VTip± oc ip±, so the 
eigenstate exhibits no net amplification nor dissipation 
(l s ±l = !)■ I n the broken-symmetry phase, ip± is not 
itself PT-symmetric but the pair satisfies VT by trans- 
forming into each other: VTip± = tp^i wnere s ± = V s ^- 
Each eigenstate in the pair spontaneously breaks VT 
symmetry; one exhibits amplification, and the other dis- 
sipation. Similar to transitions in Hamiltonian systems, 
such as the transition from real to complex eigenvalues 
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FIG. 1: (Color online) Semilog plot of S-matrix eigenvalue 
intensities log 10 |s±| 2 versus frequency uj (solid curves), for a 
ID PT-symmetric system of length L with balanced refractive 
index n = 3 ± 0.005i in each half (inset). The VT symmetry 
is spontaneously broken at ui c ~ 1418. 21/L. Dashed curves 
show the minimum scattered intensity for equal-intensity in- 
put beams of variable relative phase; in the broken-symmetry 
regime, net amplification always occurs. 

of PT-symmetric Hamiltonians, the scattering transition 
can be induced by tuning the parameters of S(lu) 12J. 

Having shown that a PT-breaking transition can take 
place in the scattering behavior of 'PT-symmetric sys- 
tems, we turn now to some concrete examples to see how 
it occurs. First, consider an arbitrary ID system with 
PT-symmetric n(x); its S-matrix is parameterized by 

where b = —irn/t € K; rR,ri are the reflection am- 
plitudes from right and left; and t is the (direction- 
independent) transmission amplitude. In general jr^l 7^ 
|rj,|, but their relative phase must be or ir. Although 
S depends on three real numbers, {Re(t), Im(t), b}, its 
eigenvalues only depend on two, \t\ and 6, for we can 
scale out the phase of t. One can show that the criterion 
for the eigenvalues of S to be unimodular is: 

\{r L -r R )/t\ = B(u,r) < 2. (10) 

For fixed Re{n(x)}, we write B as a function of uj and a 
T-breaking parameter r. On varying uj and/or r, violat- 
ing ( |10[ ) brings us into the broken-symmetry phase. 

Fig. [I] shows how the transition occurs by varying uj, 
in a simple slab of total length L with fixed n = hq ± ir 
in each half. The critical frequency can be shown to be 

uj c « ln(2n /r) c/tL. (11) 

The resulting "phase diagram" is shown in Fig.[2ja). The 
discrete points in the broken-symmetry phase correspond 
to the CPA-laser solutions which we will soon discuss; 
these lie along a line given by the equation [7] 

u « uj c + (c/tL) ln[(n 2 + l)/(ng - 1)]. (12) 



3 



(a) 



5.1 



x 10 




CPA-laser solutions 
/ 



_ 5 
c 

I 4.9 
II 

H 4.8 
4.7 

1440 1450 1460 1470 1480 1490 1500 

(b) _ 3 " L 

T xlO 



Symmetry-broken phase 



Symmetric phase 



„ 6.5 
c 

E 6 
II 5.5 



4.5 



l 

Neutral 
Loss Gain 



1440 1450 1460 1470 1480 1490 1500 
COL 



FIG. 2: (Color online) Phase diagram of S-matrix eigenvalues 
for ID 'PT-symmetric systems. The structure in (a) is the 
same as in Fig. [I] The structure in (b) has a real-index region 
of length L/10 in the center, with all other parameters the 
same as in (a). It shows re-entrant behavior as cj is varied. 
White areas show the PT-symmetric phase, and grey areas 
show the PT-broken phase; circles in (a) show the CPA-laser 
solutions. In (a), the dashed line shows the approximate phase 
boundary given by Eq. the points A, B, and C on this 

boundary are referenced in Fig. [3] For fixed ui, stronger T- 
breaking always favors broken symmetry. 



The strong oscillations of log \ s±\ in Fig. [I] occur as the 
system crosses this curve. Without fine-tuning n , the 
system will not hit one of the CPA-laser solutions exactly; 
hence, varying lu does not bring log \s±\ to ±00, though 
it can become quite large. 

With different n{x) 1 more complicated phase diagrams 
are possible. Fig. shows the effect of a real-index 
region inserted between the gain and loss regions. As we 
increase w, internal resonances in the real- index region 
cause the trapping time to oscillate with uj; hence the 
effective T-breaking perturbation strength oscillates, and 
the system re-enters the symmetric phase periodically. 

In these ID structures, the symmetry-breaking transi- 
tion can be probed by measuring the total intensity scat- 
tered by equal-intensity input beams of variable relative 
phase 4>. In the symmetric phase, two particular values of 
4>(uj, n) correspond to unimodular S-matrix eigenstates; 
for other values of <f>, both net amplification and dissipa- 
tion can be observed in the total output IS 1 ^ 2 , since the 
system is not unitary. In the broken-symmetry phase, 
where \s+\ > 1, it can be shown that 
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2) > 1, 
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where a = i(>f'ip/' t l J T' t f ; i- Thus all values of 4> for balanced 
inputs ip give rise to amplification (Fig. [I]). This signa- 
ture of the VT transition should be easily observable. 

Fig. [3] shows the behavior of the poles and zeros in the 
complex plane for the simple 2-layer structure. When 
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FIG. 3: (Color online) (a) Trajectories of S-matrix poles and 
zeros for the two-layer system of Fig. [I] Filled circles and 
squares show the zeros and poles for r = Im[n] = 0; solid 
curves show their trajectories as r is increased. Filled stars 
show CPA-laser solutions where zeros and poles meet on the 
real axis. For the middle set of trajectories, this happens at 
ujL = 1468.7, t = 4.982 x 1(T 3 . The transition occurs at the 
"anti-crossings", indicated by the representative points A, B, 
and C, which match the values of {Re[w], r} marked in Fig. [2] 
(b) Log of the transmittance for r = 0.005 and r = 0, showing 
the doubled free spectral range in the CPA-laser. 



t = 0, the poles and zeros of the T-symmetric S-matrix 
are located symmetrically around the real axis, and can 
be labeled by the parity of the corresponding eigen- 
vectors. As r increases from zero, VT symmetry re- 
quires that the zeros and poles move symmetrically in 
the complex plane. Neighboring zeros (poles) approach 
each other pairwise and undergo an anti-crossing, caus- 
ing strong parity mixing. The horizontal motion of the 
poles and zeros corresponds to the symmetric phase of 
the scatterer, and the vertical motion to the broken- 
symmetry phase. The anti-crossing corresponds to the 
phase boundary, as can be seen by comparing the points 
labeled A, B and C in Figs. [2] and [3] 

CPA-Laser - The CPA-laser points occur when a pole 
and a zero of the S-matrix coincide on the real axis, as 
shown in Fig. [3] This corresponds to the singular case 
in which \s n \ — > 0, implying |l/s*| — >■ oo, while their 
product remains unity. Physically, the scattering system 
behaves simultaneously as a laser at threshold and a CPA 
9J . Fig. [3] also indicates a very interesting property of 
these solutions: only half the zero-pole pairs flow upwards 
and reach the real axis, while the other half go off to in- 
finity and do not produce laser-absorber modes; thus, the 
CPA-laser lines have twice the free spectral range of the 
passive cavity resonances. For ujL 1, they occur at 
frequencies ui m ~ (2m + 0.5)7r/(noL), where m is an in- 
teger, exactly half-way between alternate pairs of passive 
cavity resonances. This property of the CPA-laser should 
be straightforward to demonstrate experimentally. 

In a general 'PT-symmetric system, the poles and zeros 
are related by complex conjugation, so any VT system 
that is a laser must be a CPA-laser. Even away from the 
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FIG. 4: (Color online) Semilog plot of S-matrix eigenvalue in- 
tensities log 10 |s| 2 versus frequency w, for a 2D PT-symmetric 
disk of radius R and refractive indices 1.5±0.1i. The S-matrix 
is truncated to angular momenta \m\ < 20. Most of the 41 
eigenvalues are unimodular at these frequencies; those with 
higher average m, near resonances of the passive disk, show 
greater symmetry breaking. The colors distinguish several 
different pairs of PT-broken eigenvalues. 

CPA-laser singularities, the PT-symmetric cavity in the 
broken symmetry phase is a unique interferometric am- 
plifier: for coherent input radiation not corresponding 
closely to the damped S-matrix eigenvector, amplifica- 
tion takes place — in particular, this typically occurs for 
one-sided illumination. However, coherent illumination 
from both sides, with appropriate relative phase and am- 
plitude corresponding to the damped eigenvector, leads 
instead to strong absorption. The CPA-laser is the ex- 
treme case, where the gain/loss contrast is infinite. 

The possibility of lasing in a PT system is not obvi- 
ous, as naively one might argue that photons traversing 
the device should experience no net gain. The presence 
of spontaneous symmetry breaking invalidates this argu- 
ment. Photons in amplifying modes spend more time 
in the gain region than the loss region; hence, as the 
gain/loss parameter r is increased, it eventually becomes 
possible for one mode to overcome the outcoupling loss 
and lase. As a self-organized oscillator, the CPA-laser 
will automatically emit in the amplifying eigenstate. 

Disk VT-scatterer - Consider any 2D PT-symmetric 
body in free space. Scattering in this system is naturally 
described using angular momentum channel functions, 

u%™\r,4>,u)=HZ(cjr)e ±im *, m = 0,±l,--- (14) 

For a region of linear dimensions R and average index 
n, only states with m < nkR are significantly scattered, 
so we can truncate the infinite S-matrix to N > nkR 
channels. As this S-matrix is tuned, multiple pairs of 
eigenstates can undergo PT-breaking, so the phases of 
the S-matrix are indexed by the number of pairs of PT- 
broken eigenstates. In Fig. |4j we show this behavior for 
a disk with balanced semicircular gain and loss regions, 
with the S-matrix calculated using the R-matrix method 
[10\ 113], The symmetry-breaking is non-monotonic in w, 
similar to the ID example of Fig{2jb). The regions of 
strong PT-breaking are roughly associated with the res- 



onances of the T-symmetric disk, as symmetry-breaking 
is enhanced by the dwell time in the medium. This points 
to the possibility of microdisk based amplifier-absorbers. 

Conclusion - We have shown that PT-symmetric op- 
tical scattering systems generically display spontaneous 
symmetry breaking, with a unimodular phase where the 
S-matrix eigenstates are norm-preserving, and a broken 
symmetry phase in which they are pairwise amplifying 
and damping with reciprocal eigenvalue moduli. This 
PT-breaking transition can be tested experimentally, us- 
ing ID hetero junction geometries (with realistic values 
of the gain/loss parameter) that are distinctly different 
from the paraxial geometries previously suggested for ob- 
serving PT symmetry breaking in optics [3HS]. In our 
analysis, we have neglected the role of the noise due to 
amplified spontaneous emission, which may be significant 
at the singular CPA-laser points 8]. However this noise 
should not preclude observation of experimental signa- 
tures of CPA-lasing, e.g. the doubling of the free spectral 
range relative to the passive cavity. Moreover one may 
study the intcfcromctric amplifying behavior in the bro- 
ken symmetry phase well below the CPA-laser points, 
where the noise will be substantially smaller. 
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